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Abstract. Applying the HSD transport approach to charmonium dynamics within 
the 'hadronic comover model' and the 'QGP melting scenario', we show that the 
suppression pattern seen at RHIC cannot be explained by the interaction with baryons, 
comoving mesons and/or by color screening mechanism. The interaction with hadrons 
in the late stages of the collision (when the energy density falls below the critical) 
gives a sizable contribution to the suppression. On the other hand, it does not account 
for the observed additional charmonium dissociation and its dependence on rapidity. 
Together with the failure of the hadron-string models to reproduce high V2 of open 
charm mesons, this suggests strong pre-hadronic interaction of cc with the medium at 
high energy densities. 
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1. J/^ production vs suppression in different theoretical scenarios 

The microscopic Hadron- String-Dynamics (HSD) transport calculations (employed here) 
provide the correct space-time geometry of the nucleus-nucleus reaction and a rather 
reliable estimate for the local energy densities achieved. The energy density e{r; t) - 
which is identified with the matrix element T°°(r; t) of the energy momentum tensor in 
the local rest frame at space-time (r, t) - reaches as high as 30 GeV/fm^ in a central 
Au+Au collision at = 200 GeV [1]. 

According to present knowledge, the charmonium production in heavy-ion 
collisions, i.e. cc pairs, occurs exclusively at the initial stage of the reaction in primary 
nucleon-nucleon collisions. The parametrizations of the total charmonium cross sections 
{i = Xc, J/^ 1 from NN collisions as a function of the invariant energy ^/s used in 
this work are taken from [21 El IH [S] . We recall that (as in Refs. [U El El El E]) the charm 
degrees of freedom in the HSD approach are treated perturbatively and that initial hard 
processes (such as cc or Drell-Yan production from A^A^ collisions) are 'precalculated' 
to achieve a scaling of the inclusive cross section with the number of projectile and 
target nucleons as x At when integrating over impact parameter. For fixed impact 
parameter h the cc yield then scales with the number of binary hard collisions Nun [cf- 
Fig. 8 in Ref. [5j). 

In the QGP 'threshold scenario', e.g the geometrical Glauber model of Blaizot et 
al. [9] as well as the percolation model of Satz [10], the QGP suppression '(i)' sets in 
rather abruptly as soon as the energy density exceeds a threshold value Ec, which is a 
free parameter. This version of the standard approach is motivated by the idea that the 
charmonium dissociation rate is drastically larger in a quark-gluon-plasma (QGP) than 
in a hadronic medium [TO] . 

On the other hand, the extra suppression of charmonia in the high density phase 
of nucleus-nucleus collisions at SPS energies [HI [121 lEl [E] has been attributed to 
inelastic comover scattering {cf. [3 El E [IS [ISl m [HI [HI [201 [2l] and Refs. therein) 
assuming that the corresponding J/^^-hadron cross sections are in the order of a few 
mb [22l [23l [2^ [25]. In these models 'comovers' are viewed not as asymptotic hadronic 
states in vacuum but rather as hadronic correlators (essentially of vector meson type) 
that might well survive at energy densities above 1 GeV/fm^. Additionally, alternative 
absorption mechanisms might play a role, such as gluon scattering on color dipole states 
as suggested in Refs. [261 EH [281 EH] or charmonium dissociation in the strong color fields 
of overlapping strings [6]. 

The explicit treatment of initial cc production by primary nucleon-nucleon collisions 
and the implementation of the comover model - involving a single matrix element Mq 
fixed by the data at SPS energies - as well as the QGP threshold scenario in HSD 
were explained in Ref. [H [15] (see Fig. 1 of Ref. [15] for the relevant cross sections). 
We recall that the 'threshold scenario' for charmonium dissociation is implemented as 
follows: whenever the local energy density e{x) is above a threshold value Ej (where the 
index j stands for J/\E', x^, \E''), the charmonium is fully dissociated to c + c. The 
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Figure 1. The J/^* nuclear modification factor Raa for Au + Au collisions at 
■\/s — 200 GeV as a function of the number of participants Npart in comparison to the 
data from [10] for midrapidity (full circles) and forward rapidity (full triangles). HSD 
results for the QGP 'threshold melting' scenarios are displayed in terms of the lower 
(green solid) lines for midrapidity J/^f's {\y\ < 0.35) and in terms of the upper (orange 
dashed) lines for forward rapidity (1.2 < y < 2.2) within different recombination 
scenarios (see text). The error bars on the theoretical results indicate the statistical 
uncertainty due to the finite number of events in the HSD calculations. Predictions 
for the ratio Bf^^{^')ami / B^^{J /'^)(7j/qi as a function of the number of participants 
Npart are shown in the lower set of plots. The figure is taken from [T]. 



default threshold energy densities adopted are Si = 16 GeV/fm^ for J/"^, 62 = 2 
GeV/fm^ for Xc and £3 = 2 GeV/fm^ for Two more scenarios were implemented 
similarly to the 'comover suppression' and the 'threshold melting' by adding the only 
additional assumption - that the comoving mesons (including the D-mesons) exist only 
at energy densities below some cut-energy density ecut, which is a free parameter. We 
set ecut = 1 GeV/fm^, i.e. of the order of critical energy density. 

2. Comparison to data 

In the following, we compare our calculations to the experimental data at the top 
RHIC energy of ^/s = 200 GeV. We recall that the experimentally measured nuclear 
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Au+Au, s'^^=200 GeV, comover absorption | ^ 




part 



Figure 2. Same as Fig. 4 for tlie 'comover absorption scenario' including the 
charmonium reformation channels without cut in the energy density (l.h.s.) and with 
a cut in the energy density ecut — 1 GeV/fm^ (see text for details). The figure is taken 
from ijj. 



modification factor Raa is given by 

dNiJ/^)AA/dy 
N,oU-dN{J/^)Jdy' 

wliere dN{J /^)aa/ dy denotes tlie final yield of J/\I/ in AA collisions, dN{J/'^)pp/dy is 
the yield in elementary pp reactions, Ncou is the number of binary collisions. 

Due to very high initial energy density reached (corresponding to T ^ 2Tc), 
in the threshold melting scenario all initially created J/\E', \E'' and Xc mesons melt. 
However, the PHENIX collaboration has found that at least 20% of J/\E' do survive at 
RHIC [30]. Thus, the importance of charmonium recreation is shown again. In HSD, 
we account for J/\E' recreation via the DD annihilation processes as explained in detail 
in [U |15] . Note that in our approach, the cross sections of charmonium recreation in 
D+D J/ "$+77163071 processes is fixed by detailed balance from the comover absorption 
cross section J/'if + 7neson —> D + D. But even after both these processes are added to 
the threshold melting mechanism, the centrality dependence of the RAAiJ/"^) cannot 
be reproduced in the 'threshold melting' scenario, especially in the peripheral collisions 
(see Fig. [1]). This holds for both possibilities: with (r.h.s. of Fig. [1]) and without (center 
of Fig. [1]) the energy density cut ecut, below which D-mesons and comover s exist and 
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can participate in D + Z) ^ J/\E' + meson reactions. 

Comover absorption scenarios give generally a correct dependence of the yield on 
the centrality. If an existence of D- mesons at energy densities above 1 GeV/fm^ is 
assumed, the amplitude of suppression of J/"^ at mid-rapidity is also well reproduced 
(see the line for 'comover without ec^t' scenario in Figj2l l.h.s.). Note that this line 
correspond to the prediction made in the HSD approach in pT]. On the other hand, 
the rapidity dependence of the comover result is wrong, both with and without ecut- If 
hadronic correlators exist only at e < ecut, comover absorption is insufficient to reproduce 
the J/\E' suppression even at mid-rapidity (see Fig. [2l r.h.s.). But its contribution to 
the charmonium suppression is, nevertheless, substantial. The difference between the 
theoretical curves marked 'comover + ecut^ and the data shows the maximum possible 
supression that can be attributed to a deconfined medium. 



3. Summary 

The formation and suppression dynamics of J/"^, Xc and \E'' mesons has been studied 
within the HSD transport approach for Au + Au reactions at ^/s = 200 GeV. Two 
currently discussed models, i.e. the 'hadronic comover absorption and reformation' 
model as well as the 'QGP threshold melting scenario' have been compared to the 
available experimental data. 

We find that both 'comover absorption' and 'threshold melting' fail severely at 
RHIC energies [1]. The failure of the 'hadronic comover absorption' model goes in line 
with its underestimation of the collective fiow f 2 of leptons from open charm decay as 
investigated in Ref. p2]. This suggests that the dynamics of c, c quarks at this energy 
are dominated by strong pre-hadronic interaction of charmonia with the medium in 
strong QGP (sQGP), which cannot be modeled by 'hadronic' scattering or described 
appropriately by color screening. 

On the other hand, the interaction of J/\l/ with hadrons in the late stages of the 
collision (when the energy density falls below the critical) gives a sizable contribution 
to its suppression. This contribution should not be neglected, when comparing possible 
models for QGP-induced charmonium suppression to the experimental data. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
\12] 



O. Linnyk, E. L. Bratkovskaya, W. Gassing, and H. Stocker, arXiv:0705.4443[nucl-th]. 

W. Gassing and E. L. Bratkovskaya, Phys. Rep. 308, 65 (1999). 

W. Gassing, E. L. Bratkovskaya, and S. Juchem, Nucl. Phys. A674, 249 (2000). 

E. L. Bratkovskaya, W. Gassing, and H. Stocker, Phys. Rev. C67, 054905 (2003). 

W. Gassing, E. L. Bratkovskaya, and A. Sibirtsev, Nucl. Phys. A691, 753 (2001). 

J. Geiss, G. Greiner, E. L. Bratkovskaya, W. Gassing, and U. Moscl, Phys. Lett. B447, 31 (1999). 

W. Gassing and E. L. Bratkovskaya, NucL Phys. A623, 570 (1997). 

W. Gassing and G. M. Ko, Phys. Lett. B396, 39 (1997). 

J. P. Blaizot and J. Y. Olhtrault, Phys. Rev. Lett. 77, 1703 (1996), hep-ph/9606289. 
H. Satz, J. Phys. G32, R25 (2006). 

NA50, M. G. Abreu et ai, Phys. Lett. B410, 337 (1997). 
NA50, M. G. Abreu et ai, Phys. Lett. B477, 28 (2000). 



Charmonium dynamics in heavy ion collisions 



6 



[13] NA50, M. C. Abrcu et at, Phys. Lett. B450, 456 (1999). 

[14] NA60, A. Foerster et ai, J. Phys. G32, S51 (2006), nucl-ex/0609039. 

[15] O. Linnyk, E. L. Bratkovskaya, W. Gassing, and H. Stocker, Nucl. Phys. A786, 183 (2007). 
[16] N. Armesto and A. Capella, Phys. Lett. B430, 23 (1998), hep-ph/9705275. 

[17] R. Vogt, Phys. Rep. 310, 197 (1999). 

[18] C. Gerschel and J. Hiifner, Ann. Rev. Nucl. Part. Sci. 49, 255 (1999). 
[19] D. E. Kahana and S. H. Kahana, Prog. Part. Nucl. Phys. 42, 269 (1999). 
[20] G. Spieles et ai, J. Phys. G25, 2351 (1999), Phys. Rev. C60 (1999) 054901. 
[21] L. Gcrland, L. Frankfurt, M. Strikman, H. Stocker, and W. Greiner, Nucl. Phys. A663, 1019 
(2000). 

[22] K. L. Haglin, Phys. Rev. C61, 031903 (2000). 

[23] Z. Lin and G. M. Ko, Phys. Rev. C62, 034903 (2000). 
[24] Z. Lin and G. M. Ko, J. Phys. G27, 617 (2001). 

[25] A. Sibirtsev, K. Tsushima, and A. W. Thomas, Phys. Rev. C63, 044906 (2001). 

[26] B. Zhang, G. M. Ko, B.-A. Li, Z. Lin, and B.-H. Sa, Phys. Rev. C62, 054905 (2000). 

[27] L. Grandchamp and R. Rapp, Phys. Lett. B523, 60 (2001), Nucl. Phys. A 709 (2002) 415. 

[28] D. Blaschke, Y. Kahnovsky, and V. Yudichev, Lect. Notes Phys. 647, 366 (2004), hep-ph/0410338. 

[29] M. Bedjidian et al, hep-ph/0311048. 

[30] PHENIX, H. Biisching et al. Nucl. Phys. A774, 103 (2006). 

[31] E. L. Bratkovskaya, A. P. Kostyuk, W. Gassing, and H. Stocker, Phys. Rev. C69, 054903 (2004). 
[32] E. L. Bratkovskaya, W. Gassing, H. Stocker, and N. Xu, Phys. Rev. C71, 044901 (2005). 



